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We illustrate that both of massless pseudoscalar meson fields and constitutent quark fields are 
dynamical field degrees of freedom in energy region between chiral symmetry spontaneously broken 
• scale and quark confinement scale. This physical configuration yields a renormalization running 

\ axial coupling constant g A (/j,). The one-loop renormalization prediction g A (fJ, = rn N ) = 0.75 and 

£vQ ■ g A (n = m n ) = 0.96 agree with f3 decay of neutron and 7r° — > 77 respectively. We also calculate 

I ' p ± — > 77r and lu — > -y-K decays up to the next to leading order of iV^ 1 expansion and including 

all order effects of vector meson momentum expansion. The results strongly support prediction of 
renormalization running, g A (fi = m p ) — 0.77. 

^ \ 14.40.-n,12.39.-x,11.40.Ha,13.20.Jf 

t-H ' The asymptotic freedom, chiral symmetry spontanously broken(CSSB) and quark confinement (QC) are essential fea- 
tures of QCD, which are related to the non-Abelian guage symmetry structure of QCD. They divide QCD-description 
■ of strong interaction into three different energy regions, ft is well-known that the dynamical field degrees of freedom 
in the energy region above CSSB scale are current quarks and gluons, and in the energy region blow QC scale are 
pseudoscalar meson fields purely. In the energy region between CSSB scale and QC scale, the CSSB implies that 
the dynamical field degrees of freedom should be constituent quarks (quasi-particle of quarks), gluons and Goldstone 
bosons (massless pseudoscalar meson fields) associated with CSSB. Nevertheless, it is still debated whether pseu- 
doscalar mesons are double counted when they are regarded as fundamental fields, as well as light bound states of 
quarks |j| . Theoretically, it is clear that there is not double counting problem, since fundamental pseudoscalar mesons 
are light bound states of light current quarks instead of heavy constituent quarks. Practically, however, due to lack of 
essential evidences to confirm this point, the pseudoscalar mesons are still regared as composited fields of constituent 
quarks in many works of low energy QCD in the literatures The purpose of this present paper is to provide 

an evidence to confirm that pseudoscalar mesons are fundamental dynamical degrees of freedom in the energy region 
below CSSB scale indeed. 

Our analysis is based on to solve a puzzle which also lies low energy strong interaction. In QCD below CSSB scale, 
the axial weak current with light flavors can be defined with a non-unit coupling constant, 



X 



i5/x = 5A 7 / , 7M75yV', (a = 1,2,3), 

where ip ~ (u, d) and r a are Pauli matrices. If we recognize that the wavefunction of nucleon has 5'(7(3) co i or x 
S'(7(2)fl avor symmetry, the beta decay of neutron requires g A = 0.75 Meanwhile, the PC AC (partially conserved 
axial current) hypothesis and Adler-Bell-Jackiw anomaly p] imply axial current divergence 

9"$ = y W ^° + ^9 A a c . m .^F, v F a0 , (1) 

where f„ = 185.2MeV is the pion decay constant and — d^A^ — d u A^ is photon field strength. Thus in chiral 
limit, the transition matrix element for the ir° — > 77 reduces to 
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< 7 (ci, ki), 7 (e 2 , k 2 )\n°(q) >= {2n)H 4 {q -k x - k 2 )(-i)—^g A a e . m e^e llx k lv e 2a k 2 p. (2) 

The experimental data for r(7r° — > 77) yields axial coupling constant g A = 1 for iV c = 3. This value is very different 
from one determined by (3 decay of neutron. Then how to resolve this puzzle is a open question. In this paper, we will 
provide a possible explanation on this problem in the framework of constituent quark model. In this interpretation, 
the psudoscalar mesons must be regarded as independent dynamical field degrees of freedom in this energy region. 
The renormalization of pseudoscalar-constituent quarks loop effects will lead to aixal coupling constant g A running 
with renormalization scale. Our results for one-loop renormalization are that taking g A (fi = mjv) = 0.75 as input, we 
predict g A (p — m^) — 0.96 which matchs with requirement of ir° — > 77 decay well. This result is an evdience that 
pseudoscalar mesons are fundamental fields in energy scale below CSSB. In order to confirm this conclusion further, 
we will calculate the anomal decays p ± — > 77r ± and v — > 777°. In this type of decays, effects of g A are the leading 
order. Our calculations will be up to the next to leading order of -/V" 1 expansion and including all order effects of 
vector meson momentum expansion. The results strongly support prediction of one-loop renormalization for axial 
coupling constant, g A (fi — m p ) = 0.77. 

The simplest version of chiral quark model was originated by Weinberg B , and developed by Manohar and Georgi 
jjj provides a QCD-inspired description on the constituent quark model. At chiral limit, it is parameterized by the 
following SU(3) V invariant chiral constituent quark lagrangian 

F 2 

C x = i${p + f + g A fas-)1> - + —TrfiV^UV^}. (3) 

Here TV/ denotes trace in SU(3) flavour space, ip = (u, d, s) are constituent quark fields. The A M and r M are defined 
as follows, 

= Itfidp - ir^ + ^ - il^}, (4) 
and covariant derivative are defined as follows 

V M C/t - d M U* - ilpU* + iU^r^ = -2$ , (5) 

where = + a M and r M = — a M are linear combinations of external vector field and axial-vector field 
a^, £ associates with non-linear realization of spontanoeusly broken global chiral symmetry G = SU(3)l x SU(3)r 
introduced by Weinberg M , 

e(*)^.9^(*)^(*) = M*)£(*)sL 9L,9 R eG, h($)eH=SU(3) v . (6) 

Explicit form of £($) is usual taken 

£(*) = exp{i\ a <S> a (x)/2}, [/($)= ^ 2 ($), (7) 

where the Goldstone boson $ a are treated as pseudoscalar meson octet. The constituent quark fields transform as 
matter fields of SU(3) V , 

4> — ► h($)ip, $ — ► 4>h\$>). (8) 

A M is SU(3) V invariant field gradients and T p transforms as field connection of SU(3) V 

a m — > ^($)A M fct($), r M — ► fc($)i>t(*) + h($)d^(<s>). (9) 

Thus the lagrangian ( ^) is invariant under Ggiobai x Gi oca i- 

Let us consider one-loop effects renormalization of lagrangian ([}]) In terms of defining the following "renormal- 
ization" quantities, 



l It should be pointed out that, this model as low energy effective model of QCD is not completely renormalizable. The reason 
is that more and more new divergent terms will apear when loop effects are included. Fortunately, the terms in lagrangian (^) 
are still renormalizable. 
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ip ~* i>n = z i> V2 V'> 9a -> 9ra = Z a9 A , 

m — > m R = m + 5m, F —> F R = Z^ 2 F, 

f r^ = z- 1 r M , (io) 



the divergences yielded by loop-effects of constituent quarks and massless Goldstone fields can be cancelled by divergent 
constants Z^, Z g , 5m, Z F and Z^. In this paper, we focus our attention on renormalization of axial constant g A . 
The diagrams in fig. 1 concern our purpose. In MS scheme, explicit calculation gives 

2m 2 2m 2 
= 1 - -jf9 2 RA Ne, Z^5m = -j^gLN., 

TV? TYl 3 5 

Z i> Z i = 1 " -^SL N « Z * Z * Z 9 l = 1 - T# (o + ^9 2 RA )N e , (11) 



where A x = 2irf n ~ 1.2GeV is CSSB scale and 



AT e = lim - - 7 P +ln47r. (12) 




FIG. 1. One-loop renormalization of lagrangian (g). The solid line is constituent quarks, 
dot line denotes massless pseudoscalar meson fields and wave line denotes external vector 
or axial-vector fields. 



Then renomalization lagrangian is written as follow 



2 m 2 ii 2 tti 4 

C R = i$ R (P + f R )^ R - {1 - ^g 2 RA In ^ + {^)}mM h 



m 2 3 4 ii 2 m 4 - F 2 

+U + Af ( 2 ~ 3 5 - } ln ^ + 0(^)}9 R J R ^ R + -^Tr } {V,U R ,V»U R }, (13) 

where fj, is renomalization scale. The /^-independence of renormalization requires the constituent quark mass m R and 
axial constant g RA should be /^-dependent, and should satisfy(for sake of convenience, we omit the subscript "R" in 
m R and g RA hereafter) 



where 



™<y) , 2m 2 _ 2 v 2 n ,m 4 

^y = 1+ A| ( 2"3^ )ln ^ + ° ( Ai ) ' (W) 



mjfj,) + m(fi') _ _ g A (/i) + g A (//) 

o j 9a — o ■ \ 10 ) 
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Then inputting m(/i = m p ) — 480MeVR and g A (n = m^) = 0.75, we have 



m(fi = m N ) = 494MeV, m(n = m n ) = 368MeV 

9 A (t 1 = m p) = °- 77 > 9a (M = m 7r) = 0.96. (16) 

In particular, <? A (/i = m^) — 0.96 is agree with ir° — > 77 requirement. 

Several remarks are necessary here. 1) The renormalization running of m and g A is rather larger. This point can 
be easily understood, since there is no a rather small parameter to suppress loop effects in this energy region. This 
is a feature of low energy QCD. For example, a low energy coupling constant of chiral perturbative theory at 0(p 4 ), 
L 5 , is equal to (2.2 ± 0.5) x 10~ 3 at energy scale fi = m v but (1.4 ± 0.5) x 10~ 3 at energy scale // = m p . 2) If 
pseudoscalar mesons are not independent dynamical field degrees of freedom, all diagrams in fig.l will be absent. 
Then we can not understand why the value axial constant g A should be equal to 0.75 in /3 decay of neutron but be 
equal to 1 in 7r° — > 77. Therefore, our results in eqs. (14|) and ( |l3| ) is a direct support that massless pseudoscalar 
mesons as Goldstone bosons associating CSSB are fundamental dynamical degrees of freedom. 3) The loop effects 
also generate some high derivative terms contributing to axial constant, e.g., ^d 2 ^^. However, we expect that 
this type terms are suppressed by p 2 / A 2 , expansion (where p denote four-momentum of axial current). For (3 decay of 
neutron and 7r° — > 77, the transition momentum are both very small comparing with K x . So the contributions from 
high derivative terms are very small that it can not explain the difference of axial constant in two reactions. 4) If 
we inserted Adler-Bell-Jackiw anomal term with unit axial constant into original lagrangian (^|), unambiguously, it 
should cause double counting, since triangle diagrams of constituent quarks also generate Adler-Bell-Jackiw anomally. 
5) In lagrangian (^) pseudoscalar mesons are massless. It is consistent with Goldstone theorem for CSSB. 

In the following, for checking the result of eq. ( |l4| ) , we will provide a complete calculation on decays — > 77r ± and 
lu — ► 77T°. Since the chiral expansion converges slowly at energy scale of vector meson masses we have to include 
high order contribution of chiral expansion into our results. In this paper, our calculation will be up to the next to 
leading order of iV c _1 expansion, and including all important effects of momentum expansion. 

Due to WCCWZ realization for vector meson resonances H[|, the lagrangian can be easily extended to included 
the lowest vector meson resonances 



C x = i$(p + f + g A ^l 5 - TO ~ + ^TrfiV^UVtf} + ^m 2 Tr f {V^}. (17) 

Here denotes vector meson octet and singlet, or more convenience, due to OZI rule, they are combined into a 
singlet "nonet" matrix 



It transforms homogeneously under SU(3) V 



V k*- k;° ^ j 



(18) 



Vn^hityVuhti*), (19) 

Thus the lagrangian ( |l7|) is still invariant under Ggiobai x Gi oca i. The effective action describing meson interaction 
can be obtained via integrating over degrees of freedom of fermions in lagrangian ( |l7j ) 

e iSett = VqDqeI d4xCx(x) =< vac,out\in,vac> v>A ,r, (20) 



where < vac, out\in, vac >y.A,r is vacuum expectation value in presence external sources. The above path integral 
can be performed explicitly, and heat kernal method Jl(],[ll| has been used to regulate the result. However, this 
method is extremely difficult to compute very high order contributions in practice. This diffculty can be overcomed 



2 The 0(p 4 ) chiral coupling constant L$ = (|g 2 — ) j§~ ls predicted in this model. Here g and Bo are two constants 
absorbing divergence from quark loops. In particular, g — ty^ 1 and Bo = m%/{m u +md) — 2GeV is fitted at /1 = m p Thus 
m = 480MeV is fitted by L 5 (p = m p ) = (1.4 ± 0.5) x 10" 3 . 
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via calculating one-loop diagrams of constituent quarks directly. This method can capture all high order contributions 
of the chiral expansion. 

In interaction picture, the equation( EG) is rewritten as follow 



. f ,4 d4 P2 d 4 p n - 4 

= J d pi^-y I ---^-y i n„(pi,---,p„) ( 5 [PI-P2 Pn) 



n=2 



where T q is time-order product of constituent quark fields, C l x is interaction part of lagrangian( |l7|) , fl n (pi, ■ ■ ■ ,p n ) is 
one-loop effects of constituent quarks with n external sources, P\,pi, • • • ,Pn are four-momentas of n external sources 
respectively and 

n n (Pl,"--,Pn-l) = J d 4 p n fi n (pi,- ■ ■ ,p n )8 4 (pi - p 2 p n )- ( 22 ) 

To get rid of all disconnected diagrams, we have 

oo 

S c ff — ^ ' S n , 

n=l 

5i=ni(0), (23) 

/" d 4 Pl d 4 p„-i , , , . 

A " = / 7^1 /r, n 4 n ra (pi, ■ ■ ■ ,Pn-l), (n>2). 



(2tt) 4 (2tt) 

Hereafter we will call S n as n-point effective action. 

At the leading order of N^ 1 expansion, the 3-point pjir and ujjtt effective action is generated by triangle diagram 
of constituent quarks 

n 3(g,fc) = 3w 2gf eg A J dx i J dx 2 {l - {x x + x 2 ){l - xi - x 2 ) 1 ^y 1 

xe^qMp^A^miq - k) + 3oj^q)A,Ak)ir a (q - k)}, (24) 

where g is an universal coupling constant which absorbs nonrenormalizable logarithmic divergence from constituent 
quark loops. In the above equation, g is from vector meson field normalization 

i 1 i 1 

In ref. @, .9 = tt" 1 has been fitted by the first KSRF sum rule ||§. 

It is rather complicate to calculate contribution from meson one- loop, i.e., the next to leading order contribution 
of N~ 4 expansion. The diagrams in fig. 2 concern to our calculation. Here it is requirement of unitarity of S'-matrix 
to sum over chain-like approximation in fig.2-b). 



.11 n-loop 
0), p .... K at p / + ... + m ,p> .... 



i r \ 7 

a) b) 



FIG. 2. Meson loop correction to p — > "/n and u — * jtv decay. The dot internal lines 
denote physical pseudoscalar meson fields (massless pion and massive K and r\. a) Tadpole 
loop, b) Chain-like approximation. 
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For obtain the tadpole correction in fig.2-a), at leading order of N' 1 expansion, vertices of p-fir — KK(r]rj) and 
lujtt — KK{rjrf) are needed if we treat pion as massless particle. In scheme of dimensional regularization, only the 
vertices generated by triangle diagram of constituent quarks yield non-zero contribution, 

SU 3 {q,k) = eg A J dx x J dx 2 {l - (a* + x 2 )(l - x x - x 2 )^y 1 e^ q a k^ VV 

Tr{^{p\Xq)\ l + 3w„( g ) + eQM9))(d{k)X j + 3u v (k) + eQA„{k))([X a , \ l ]\ b + X a [X l , X b }) 
+ e A„(q)([X a , Q]X b + X a [Q, A 6 ])^* (fc)A' + 3^(fc))V 

+e{p i M (q)X i + 3u M (q))A v (k)([\ a , Q]X b + X a [Q, X b })7r 1 }, I = 1,2, 3; a, b = 4, 8) (25) 

where Q = diag{2/3, —1/3, —1/3} is charge operator of quark fields, ip a denotes K— or 77-meson fields in internal line. 
For sake of convenience, we can assume that the masses of K and r\ are degenerated. Then integral over ip fields, we 
have tadpole correction to pjir and uj^/tt vertices 



nr d (9, k) = --cn 3 (q, k), c = ^W(^) e/2 r(i - y) - ^%x, (26) 

3 8w z f4 rnt 2 8tt^/^ 



where n 3 (<7, k) is defined in eq.( ^4j), A absorb the quadratic divengence from meson loops. A ~ 2/3 has been fitted 
by Zweig rule M. 



(O.P y (O.P 



\ 



+ ... 



a) b) 



FIG. 3. The chain-like approximation in fig.2-b) is equivalent to contribution of fig.3-a). 
It is generated through contracting an effective VPP(where "V" denotes p or u> meson, 
"P" denotes pseudoscalar meson) vertex and 771--PP vertex. Here the effective VPP vertex 
includes not only tree level vertex but also meson loop contribution (see fig.3-b). 



Next we will calculate the contribution in fig.2-b). It is equivalent to calculate the diagram in fig.3-a). In other 
words, we need to an effective VPP (where "V" denotes p or uj meson, "P" denotes pseudoscalar meson octet) vertex 
and 77T-PP vertex. In fig.3-b) we have shown that the effective VPP vertex includes not only tree level vertex but 
also meson loop contribution. 

At leading order of -/V" 1 expansion, the anomal 77T-PP vertex is generated by both of triangle and box diagrams of 
constituent quarks, 

r 7 ^ 3 p(<72, h,k 2 , k 3 ) = ^eB(< Z 2 )e^A AI (< Z2 )fc ll ,fc 2Q fc 3/ 3Tr / {gP(fc 1 )P(fc 2 )P(fc 3 )}, 



AT f 1 fl-xi „2 „2 

B (l ) = T1TT9a( / dx x / dx 2 {l-^(x 1 +x 2 )(l-x 1 -x 2 )}- 1 -^), (27) 

"71" Jn Jq Jq iff! 

where q 2 = {k\ + k 2 ) 2 . The effective p-PP vertex and u> — KK vertex(here we omit those vertices suppressed by 
isospin conversation) have been derived in ref. RJ 

T pP p(q, k x ,k 2 ) = \g pP p{q 2 ){q 2 k 2p - q^q ■ k 2 )Tr f {p^q)P(k 1 )P(k 2 )}, 

r uKK (q, h,k 2 ) = ^KK{q 2 ){q 2 k 2 ^ - q^q ■ k 2 )co^{(K+ (k^R- (k 2 ) + K° (fc x ) K° (k 2 )) + (h «- k 2 )}, (28) 

with 

A^q 2 ) + g 2 A A 2 (q 2 ) 



^K K (q 2 ) = „ m ; l ":; / ; A •^ ' , (29) 



s /2(l + 2C)(l + E(g a ))' 

A.jq 2 ) + g 2 A A 2 (q 2 ) 
g/2(l + 20(l-2E Jf (g2))' 



G 



where 



Ml)=9 dx-x(l-x)ln(l- -i 



7T 

l 2 (? 2 ) = -.9 2 



—| / dxi / dx 2 (xi - x 2 x 2 ){l + 
2tt 2 Jo Jo 



m 2 — Xi(l — Xi)(l — x 2 )q 2 



M^- Xl{l ~ Xl){ l~ X2)q2 )}, 



S (<?) = U + 2/2(1 + }( 4E ^(9 ) - 2S M<7 )) 5 (30) 

^(9 2 ) = ^2 + jf * • *(1 - «) In + ^(-1)% 2 - iml)}, 

Wl 2 ) = T7T272 i A K - T) + T d *K - *(! - ^ ln (1 - ^-^)}- 

Then due to soft-pion theorem, the pseudoscalar meson loops in fig.2-b) contribution to pjn and lujtt vertices as 
follow 

Ill- l ° op ( q ,k) = -eq 2 g pPP (q 2 )[A(q 2 ) + 2A(Q)][^(q 2 ) - \^ K (q 2 )]e^ q^pp^A^k^q - k) 

+ \eq 2 9ujKK (q 2 )[A{ q 2 ) + 2A{0)]Y. K (q 2 )e^ q a k fj uj p {q)A„{k)^ {q - k). (31) 

Here pr/ir coupling receives contributions from both of pion-loop and i"T-loop, but lj^tt coupling receives dominant 
contribution from JT-loop only. 

Eq.([3l|) tegother with eqs.(p4|) and ( p6| ) give "complete" p^fir and lujtt coupling up to the next to leading of iV^ 1 
expansion at least, 

ng(g, k) = ee^qakpig^iftplMMQniq - *) + g^{q 2 )u:,{q)A v (k)^ \q ~ k)}, 

g P7 7r(g 2 ) = „ 2 % g A (l ~ oO / / dx 2 {l- (Xl + X 2 )(1 -Xi - X 2 ) — — j" 1 

'J 7 '' g/lT O Jo JO ^ m 

-g 2 3pPp(g 2 )[A(q 2 ) + 2A(0)][^(q 2 ) - ^(<? 2 )], (32) 



2 "J g A (l ~ qC) / ^1/ ^ 2 {1- (Xl +X 2 )(l-Xl -X 2 )— j}" 1 

7r 2 g/7r 3 J Jo 2m 2 - 

1 g^^/f/f ((Z 2 )[A(g 2 ) + 2A(0)]Z K (q 2 ). (33) 



2' 

For <? A (/i = m p ) = 0.77, the above results yield 

B(p ± -> 71^7) = 4.88 x 10~ 4 , B(u> -► 7r°7) = 8.9%. (34) 

These results agree with data |J| B(p ± -> 71^7) = (4.5 ± 0.5) x 10~ 4 and B(w -> 7r 7) = (8.5 ± 0.5)% well. If we 
take g A = 1 as a comparison, the theoretical predictions are 

B{p ± -> 71^7) = 8.33 x 10~ 4 , B(lu -> 7T 7 ) = 15.0%. (35) 

These results obviously disagree with data. Therefore, the theoretical prediction (|34|) provides an evidence to confirm 
the result in eq. (|l4|). 

To conclude, we discuss the one-loop renormalization in Manohar-Georgi model. This renormalization leads to 
running of axial coupling constant and constituent quark mass. Inputting g A (p = mjsr) = 0.75 which agree with j3 
decay of neutron, the prediction g A (p — m n ) — 0.96 agree with reqiurement of ir° — > 77 decay. We also notice that, 
if massless pseudoscalar meson fields were not independent degrees of freedom in energy between CSSB scale and 
QC scale, we could not interpret /3 decay of neutron and ir° — > 77 decay simultaneously. For example, because axial 
constant is not unit in Extend Nambu-Jona-Lasinio model [ 
decay. The similar problem also exists in the models of ref. 



, the model can not yield right prediction for n — > 77 
y. Therefore, our result in this paper can be reagrd as 
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an evidence that pseudoscalar meson fields are independent degrees of freedom in energy between CSSB scale and 
QC scale. 

The renormalization of axial constant also predicts g A {^ = m p ) = 0.77. For confirming this prediction, we calculate 
anomal vector meson decays, p — > jtt and u> — ► 77r. In this type of decays, effects of g A is leading order. Our 
results are up to the next to leading order of iV" 1 expansion and include all order effects of vector meson momentum 
expansion. The theoretical predictions for branch ratios of these two decays agree with results of renormalization and 
is against g A = 1. 
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